Interplay of atomic displacements in the quantum magnet (CuCl)LaNb207 
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We report on the crystal structure of the quantum magnet (CuCl)LaNb207 that was con- 
troversially described with respect to its structural organization and magnetic behavior. Using 
high-resolution synchrotron powder x-ray diffraction, electron diffraction, transmission electron mi- 
croscopy, and band structure calculations, we solve the room-temperature structure of this com- 
pound [a-(CuCl)LaNb207] and find two high-temperature polymorphs. The 7-(CuCl)LaNb207 
phase, stable above 640 K, is tetragonal with Osub = 3.889 A, Csub = 11.738 A, and the space group 
PA/mmm. In the 7-(CuCl)LaNb207 structure, the Cu and CI atoms are randomly displaced from 
the special positions along the {100} directions. The /3-phase (osub x 2asub x Caub, space group 
Pbmm) and the a-phase (2asub x 2asub x Csub, space group Pham) are stable between 640 K and 
500 K and below 500 K, respectively. The structural changes at 500 and 640 K are identified as 
order-disorder phase transitions. The displacement of the CI atoms is frozen upon the 7 — >■ /3 trans- 
formation, while a cooperative tilting of the NbOe octahedra in the a-phase further eliminates the 
disorder of the Cu atoms. The low-temperature a-(CuCl)LaNb207 structure thus combines the 
two types of the atomic displacements that interfere due to the bonding between the Cu atoms 
and the apical oxygens of the NbOe octahedra. The precise structural information resolves the 
controversy between the previous computation-based models and provides the long-sought input for 
understanding (CuCl)LaNb207 and related compounds with unusual magnetic properties. 

PACS numbers: 61.66.Fn, 61.05.cp, 61.05.jm, 71.15.Mb 
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I. INTRODUCTION 

Superexchange interactions underlie the magnetic phe- 
nomena in most transition metal compounds. Such inter- 
actions involve orbitals of non-magnetic cations and are 
therefore highly sensitive to subtle changes in the crystal 
structure. This high sensitivity can be used to effectively 
vary the properties of the system by applying external 
pressure^' or by performing chemical substitutionsl^I On 
the other hand, it underlies the crucial importance of 
a precise structure determination. Inaccurate structural 
information can lead to irrelevant or even wrong conclu- 
sions regarding the magnetic properties of the system.EEl 

(CuCl)LaNb207 is a parent compound for the fam- 
ily of layered materi als, derived from perovskite-type 
Dion-Jacobson phases.^^ In such compounds, the two- 
dimensional (2D) perovskite blocks alternate with (MX) 
layers, where M^+ is a transition metal and X is CI 
or Br (Fig. [l]). The Cu-containing compounds appear 
to be especially interesting with respect to the unusual 
magnetic properties because of their low, spin-i mag- 
netic momen t th at leads to quantum behavior at low 
temperatures.'^ (CuCl)LaNb207 shows a spin gap and 
Bose- Einstein condensation of magnons in high magnetic 
fields.^^ The bromine analog (CuBr)LaNb207 under- 
goes stripe antiferromagnetic (AFM) ordering,^^ while a 
related (CuBr)Sr2Nb30io compound shows a magnetiza- 
tion plateau at one-third of the saturation. ^^^ All these in- 
teresting quantum effects remain poorly understood from 
the theoretical side due to a lack of adequate spin models 
that can be only developed, based on accurate structural 
information. 

The original study proposed (CuCl)LaNb207 to have 



tetragonal symmetry with both Cu and CI atoms on the 
four-fold axes.^ However, the high atomic displacement 
parameter of CI {Uiso = 0.13 A^) evidenced a shift of 
the CI atom from this high-symmetry position. If the 
overall tetragonal symmetry is preserved, the CI atoms 
have to be disordered over several equivalent sites. On 
the other hand, an ordered arrangement of the CI atoms 
would imply the formation of a superstructure with re- 
duced symmetry and, presumably, with an extended unit 
cell. Since the orbitals of the CI atoms mediate superex- 
change interactions, the precise determination of the CI 
position is crucial to understand the magnetic properties. 

Later studies included neutron powder diffractiorP^l 
and synchrotron x-ray powder diffraction experimentalSl 
but neither of both showed any superstructure reflections 
that could arise from an ordered arrangement of the dis- 
placed CI atoms. The superstructure reflections were ob- 
served by electron diffraction ,U^ yet no quantitative struc- 
tural information has been extracted from these data. 
Nuclear magnetic/quadrupole resonance (NMR/NQR) 
studies evidenced unique positions of the Cu, CI, and 
La atoms and confirmed the formation of an ordered 
structure.!^ The subsequent analysis of electric field gra- 
dients (EFGs) showed sizable violation of the tetragonal 
symmetry. Based on these results, Yoshida et al^^ em- 
pirically established a structural model. In their model, 
the neighboring Cu atoms are linked into a dimer via 
two CI atoms with shorter Cu-Cl bonds. The neigh- 
boring dimers are connected by longer Cu-Cl bonds only 
(the top right panel of Fig.fTl). The assignment was solely 
based on the type of the superstructure and on the hyper- 
fine field at the CI site, while no quantitative evaluation 
of EFGs or exchange couplings was performed. 
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FIG. 1. (Color online) Left panel: crystal structure of 
(CuCl)LaNb207. The Cu atoms are bonded to the CI atoms 
and to the 0(3) atoms of the NbOg octahedra. Right panel: 
the ordering models for the (CuCl) layers, as proposed in 
Refs. [T51 (top) and [IT] (bottom). Thick lines show short Cu- 
Cl bonds. The unit cell in the left panel denotes the tetragonal 
sub cell, according to Ref. [5] 



compounds. The choice of the unit cell and the crys- 
tal symmetry can also be a subtle issue. Therefore, the 
experimental input is still desirable to derive an unam- 
biguous structural model. In the following, we show 
that high-resolution synchrotron powder x-ray diffraction 
(SXRD), combined with electron microscopy and assisted 
by band structure calculations, leads to an unambiguous 
structure solution and resolves the existing controversy. 
Moreover, the variable-temperature SXRD experiments 
reveal three polymorphs of (CuCl)LaNb2 07. These poly- 
morphs are further denoted as a, /3, and 7 according to 
the temperature ranges of their stability. 

The paper is organized as follows. In Sec. |ll) we re- 
view the experimental and computational methods. In 
Sec. |III[ we determine the room-temperature (RT) a- 
(CuCl)LaNb207 structure and investigate the temper- 
ature evolution of the crystal structure. Then, we cor- 
roborate our conclusions by electron diffraction, electron 
microscopy, ab initio structure relaxation, and EFG cal- 
culation. In Sec. |IV[ we perform the symmetry analysis, 
discuss the interplay between the structures of different 
polymorphs, and the implications of our findings for the 
magnetic properties of (CuCl)LaNb2 07. 



II. METHODS 



In an earlier study, we attempted an ab initio com- 
putational approach to the problem.E^ Employing den- 
sity functional theory (DFT) calculations, we relaxed 
the crystal structure of (CuCl)LaNb207 and found that 
the system tends to undergo a Jahn- Teller-type distor- 
tion. The displacements of the CI atoms lead to two 
short and two long Cu-Cl bonds. Additionally, each 
Cu atom has two short Cu-0 bonds to the oxygens 
of the NbOg octahedra. Then, the short Cu-Cl bonds 
align along one of the intralayer directions, and chains of 
corner-sharing CUO2CI2 plaquettes are formed (bottom 
right panel of Fig. IlJ. This DFT-based model readily 
reproduced the experimental EFGs for the Cu and CI 
sites. It also explained the low energy scale of the ex- 
change couplings and the non-trivial exchange coupling 
between fourth-neighbors, evidenced by the inelastic neu- 
tron scattering.EES Still, this computation-based model 
did not yield the full solution of the problem, because it 
did not explain the large experimental asymmetry of the 
EFG at the La site. Moreover, the subtle balance of ferro- 
magnetic and antiferromagnetic superexchange couplings 
was evidenced by experimental studies,^ but could not be 
fully reproduced in band structure calculations. Further 
on, a subsequent computational studj- challenged our 
model and claimed a different structure that is similar to 
the earlier proposal by Yoshida et al^ 

Band structure calculations provide an effective ab ini- 
tio route of structure determination. However, they ex- 
perience conceptual difficulties in the treatment of elec- 
tronic correlations, relevant for many transition metal 



Powder samples of (CuCl)LaNb207 were obtained via 
a two-step procedure. First, we prepared RbLaNb207 by 
firing a mixture of La203 and Nb205 with a 25%-excess 
of Rb2C03 in air at 1050 °C for 24 h. The resulting 
white powder was washed with water to obtain single- 
phase RbLaNb207. At the second step, RbLaNb207 was 
mixed with a two- fold excess of anyhydrous CUCI2. To 
avoid any traces of moisture, the commercial anhydrous 
CuCl2 was dried overnight at 80 °C and further han- 
dled in a glove box. The mixture of RbLaNb207 and 
CuCl2 was pressed into a pellet, sealed into an evacuated 
quartz tube, and fired at 400 °C for 48 h. The resulting 
samples were washed with water to remove CuCl2 and 
RbCl. Finally, the samples were dried overnight at 80 °C. 
Note that we used a higher annealing temperature com- 
pared to the previous studies: 400 °C (673 K) vs. 325 °C 
(598 K) in Refs. |5] and [9l The higher annealing tempera- 
ture does not lead to a decomposition of (CuCl)LaNb207 
and rather improves the sample quality, because the high- 
temperature 7-modification of (CuCl)LaNb207 is formed 
during the annealing and transformed into the RT a- 
modification upon cooling (see Sec. IIIBl. Our samples 
revealed sharp superstructure refiections at RT, while 
the samples prepared at 325 °C showed a more diffuse 
scattering at certain positions. According to XRD, the 
samples were single-phase and fully reproducible. Sev- 
eral samples were additionally checked with magnetiza- 
tion measurements (MPMS SQUID, 2 - 380 K tempera- 
ture range, 0.1 — 5.0 T field range) and showed essentially 
identical susceptibility curves that were also in agreement 
with previous reports.^ 



To test the sample purity, we collected laboratory pow- 
der x-ray diffraction data using a Huber G670f Guinier 
camera (CuKd radiation, Image Plate detector). The 
high-resolution SXRD patterns were measured at the 
ID31 bcamlinc of the European Synchtrotron Radiation 
Facility (ESRF) with a constant wavelength of about 
0.4 A. The data were collected by eight scintillation de- 
tectors, each preceded by a Si (111) analyzer crystal, in 
the angle range 29 = 1 — 40 deg. The powder sample 
was contained in a thin- walled borosilicate glass capillary 
with an external diameter of 0.5 mm. To achieve good 
statistics and to avoid the effects of the preferred orienta- 
tion, the capillary was spun during the experiment. The 
sample was cooled below RT in a He-flow cryostat (tem- 
perature range 100 — 300 K) or heated above RT with a 
hot-air blower (300 - 660 K). 

The structure refinement was performed in the 
JANA2000 program. ^^ The conventional profile matching 
and the Le Bail fitting procedure were strongly impeded 
by the weak reflection splitting due to a small orthorhom- 
bic distortion. Therefore, we also used the Rietveld re- 
finement to extract the unit cell parameters at different 
temperatures. To achieve maximum accuracy, the zero 
shift and the asymmetry parameter were fixed according 
to the refinement of the RT data. The symmetry analysis 
was done in the ISDDISPLACE programPl 

The samples for transmission electron microscopy 
(TEM) were prepared by crushing the powder sample in 
ethanol and depositing it on a holey carbon grid. Selected 
area electron diffraction (ED) patterns were recorded 
using a Philips CM20 microscope. High-angle annular 
dark-field scanning TEM (HAADF STEM) images were 
taken with a Tecnai G2 microscope. High-resolution 
TEM (HRTEM) images were recorded on JEOL 4000EX 
and Tecnai G2 microscopes. The HRTEM images were 
simulated by means of the JEMS software. 

To confirm the refined crystal structure, we performed 
ab initio geometry relaxation and evaluated the EFGs 
for the a-modification, the only ordered polymorph of 
(CuCl)LaNb207. The band structure calculations were 
done with the FPLO9.00-33 codepl that performs full- 
potential DFT calculations applying the basis set of 
atomic-like local orbitals. Both the local density approx- 
imation (LDApZl and the generalized gradient approxi- 
mation (GGApS for the exchange-correlation potential 
were used. To account for correlation effects, we used 
the mean-field LSDA-|-C//GGA-|-t/ approaches with the 
around-mean-field double-counting correction schemcl^S 
The on-site Coulomb repulsion and the on-site exchange 
parameters were fixed at Ud = 5.5 eV and J^ = 1 eV, ac- 
cording to the previous studies.lJ^ The variation in Ud in 
the range of 5.5±1 eV and the type of the magnetic order- 
ing had little inffuence on the resulting geometry and on 
the EFG values. Residual forces in the relaxed structures 
did not exceed 0.01 eV/A. All the calculations were per- 
formed in the 48-atom unit cell with lattice parameters 
derived from the experimental structure refinement. The 
fc-mesh included 192 points in the symmetry-irreducible 



part of the first Brillouin zone. The general information 
on the band structure of (CuCl)LaNb207 has been re- 
ported in a previous publication.^^ 



III. RESULTS 

A. Room-temperature 
a-(CuCl)LaNb207 crystal structure 

Contrary to the previous reports,'^'^ the laboratory 
XRD powder pattern shows several weak reflections, 
violating the earlier proposed tetragonal unit cell of 
(CuCl)LaNb207 with Osub = ''sub — 3.87 A and Csub — 
11.8 A. The high- resolution SXRD experiment reveals 23 
additional reflections of this type. The resulting pattern 
can be indexed in a supercell with a ~ 6 ~ 2asub = 7.74 A 
and c ~ 11.8 A. The temperature evolution of the super- 
lattice reflections (see Sec. HIB and Fig. Isl) rules out 
their possible extrinsic origin and assigns them to the 
(CuCl)LaNb207 phase. 

The room-temperature SXRD data are consistent with 
the tetragonal symmetry and do not demonstrate any ap- 
parent reflection splitting that could originate from the 
symmetry reduction. The ED patterns, however, clearly 
indicate an orthorhombic symmetry (see Sec. |IIIC |. 
Moreover, the temperature evolution of the SXRD pat- 
tern also evidences a subtle difference between the a and 
b parameters (see Fig. [4]) and suggests an orthorhom- 
bic symmetry. At RT, the difference between a and b 
is, by chance, negligible. Therefore, no reflection split- 
tings could be observed within the present experimental 
resolution.^ As one takes into account the actual or- 
thorhombic symmetry, the space group Pham or its sub- 
group Pba2 are identified, based on the reflection condi- 
tions k — 2n for Okl and h = 2n for hOl. These reflec- 
tion conditions were also confirmed by the ED study (see 
Sec. IHC). In the following, we adopt the Pbam space 



group for the structure refinement. We did not find any 
signature of the disorder that could point to an acentric 
crystal structure. 

In (CuCl)LaNb2 07, a superstructure can be caused 
by cooperative tilts of the NbOg octahedra and/or by 
atomic displacements in the (CuCl) layers. The possible 
modes of the tilting distortion and the atomic displace- 
ments are classified by considering a virtual transforma- 
tion of the parent Osub x Osub x Csub tetragonal structure 
with the PA/mmin space symmetry to the orthorhom- 
bic 2asub X 26sub X Csub superstructure with the Pbam 
space symmetry. To denote the tilting patterns, we use 
the Glazer's notatiorPS in the form a^b'^c^, where a, 6, 
and c stand for different tilt angles around the respective 
crystal axes, +/— denote the in-phase/out-of-phase ro- 
tations, respectively, and indicates the absence of the 
tilt. The atomic displacements within the (CuCl) lay- 
ers are confined to the ab plane due to the mirror plane 
symmetry of the Pbam space group. If all the atoms 
with a certain y coordinate have the same shifts along a 
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FIG. 2. (Color online) Three models of the z = 1/2 (CuCl) layer based on the refinement of the RT SXRD data. Gray (thick) 
and yellow (thin) lines show the Cu-Cl bonds in the ranges 2.3 — 2.5 A and 2.5 — 2.8 A, respectively. The numbers under the 
plots list the refinement residuals Ri/Rp. 



(or, similarly, all the atoms with a certain x coordinate 
have the same shifts along b), the resulting displacement 
is called "in-phase". In contrast, the "out-of-phase" dis- 
placement implies the shifts of these atoms in opposite 
directions. 

Since different combinations of the atomic displace- 
ments and octahedral tilts are possible, we construct sev- 
eral starting models for the refinement. These models 
are conveniently classified by the position of the inver- 
sion center that determines the arrangement of the a and 
b glide planes, while these planes further confine possi- 
ble distortions. If an atom lies on the glide plane, the 
displacements along the glide direction should be "in- 
phase". In contrast, the glide plane passing between 
the two atoms imposes the "out-of-phase" displacements 
along the glide direction. The a and b glide planes do 
not allow for the in-phase octahedral tilts around a and 
b. The out-of-phase tilt around a is only possible if the 
a glide plane passes through the Nb atom, while the b 
glide plane contains the common oxygen atom of the two 
octahedra. Any other arrangement of the glide planes 
would completely forbid the tilts around a and b (a" 6°). 
The following starting models are possible: 

I) The inversion center is at the CI site, hence no dis- 
placements of the CI atoms are possible. The Cu and 
Nb atoms lie on the glide planes. Then, the Cu atoms 
can shift "in-phase" along both a and b axes, but the 
octahedral tilts around a and b are forbidden (the tilting 
pattern a'^b'^c'^). 

II) The inversion center is placed between the two Cu 
atoms. Then, the a glide plane contains the Cu and 
Nb atoms and passes between the CI atoms, while the b 
glide plane contains the CI atoms and passes between the 
Cu atoms. The opposite arrangement (Cu and Nb on b 
planes, CI on a planes) is also possible and leads to two 
different options: 

Ila) Left panel of Fig. [2} the CI atoms follow the "in- 
phase" displacement along b and the "out-of-phase" dis- 
placement along a. The Cu atoms show the "in-phase" 



displacement along a and the "out-of-phase" displace- 
ment along b. The tilt system is a''b''^c'^. 

lib) Middle panel of Fig. [2} the CI atoms follow the 
"in-phase" displacement along a and the "out-of-phase" 
displacement along b. The Cu atoms show the "in-phase" 
displacement along b and the "out-of-phase" displace- 
ment along a. The tilt system is aPb^c^ . 

Ill) Right panel of Fig. [2} the inversion center is at 
the Cu site, hence the displacements of the Cu atoms 
are forbidden. The CI atoms lie on the glide planes can 
be displaced "in-phase" along both a and b axes. No 
octahedral tilting is possible {aPb^c^). 

Model I can be rejected immediately: it does not al- 
low to displace the CI atoms, while the previous studies 
evidenced the pronounced shift of CI (see Sec. |l|. The 
model III implies two inequivalent Cu positions and for- 
mally violates the NMR/NQR data that suggest a single 
Cu position.liSl However, the positions with similar lo- 
cal environment are not necessarily resolved by nuclear 
resonance techniques. Therefore, we also considered the 
model III in the structure refinement. 

The models Ila and lib show a cooperative tilt of the 
NbOg octahedra around either the a or 6 axis. In the 
model Ha, the structure of the (CuCl) layer resembles the 
model by Yoshida et al^^ with pairs of Cu atoms con- 
nected into dimers via shorter Cu-Cl bonds (left panel 
of Fig. [2]). The model lib can be derived from our 
computation-based modeP^ and shows chains of corner- 
sharing CUO2CI2 plaquettes (middle panel of Fig. [2]). In 
the model Ila, the local environment of Cu is fairly irreg- 
ular with one very short (2.26 A) and one longer (2.75 A) 
Cu-Cl bond. These bonds have a cis-configuration with 
a Cl-Cu-Cl angle of 76°. In contrast, the model lib 
yields a typical CUO2CI2 plaquette with two Cu-Cl dis- 
tances of about 2.4 A in a iraras-configuration (a Cl-Cu- 
Cl angle of 173°). In the model III, the two inequivalent 
Cu atoms show different local environments, with one of 
these atoms being clearly underbonded (all the Cu-Cl 
distances are above 2.5 A). 
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FIG. 3. (Color online) Experimental, calculated, and differ- 
ence SXRD patterns for Q-(CuCl)LaNb207. Ticks show the 
reflection positions. 



The refinement residuals for the three models are sum- 
marized in Fig. [2] Based on this comparison, we select 
model lib as the correct structure of a-(CuCl)LaNb207. 
This choice is further supported by the results of band 



structure calculations (Sec. HID) and by the symmetry 
analysis (Sec. IV). It is also consistent with empirical 



crystal chemistry arguments that suggest the CUO2CI2 
plaquettes as a typical local environment of Cu^^ in 
oxychlorides .123 Atomic coordinates for the model lib are 
listed in Table llj while representative interatomic dis- 
tances are given in Table IV Figure [3] shows experimen- 
tal, calculated, and difference SXRD patterns. 

In the structure of a-(CuCl)LaNb207, the NbOe oc- 
tahedra are tilted around the b axis, while the chains of 
the corner-sharing CUO2CI2 plaquettes also align along 
the b direction. The resulting crystal structure is very 
similar to our previous DFT-based model that, for the 
sake of simplicity, did not consider the tilts of the NbOg 
octahedra.^^ The tilting distortion doubles the unit cell 
along a and is coupled to the shifts of the Cu atoms. Yet 
the local environment of Cu and the topology of the layer 
are retained. 



B. Temperature evolution of the crystal structure 

The unusual magnetic properties of (CuCl)LaNb2 07 
are observed at low temperatures, which makes the low- 
temperature crystal structure most relevant for magnetic 
studies. On the other hand, the evolution of the super- 
structure reflections above RT can give important infor- 
mation about their origin. Therefore, we investigated 
(CuCl)LaNb207 both above and below RT. 

Figure |4] shows the temperature dependence of the 
lattice parameters and the unit cell volume. The or- 
thorhombic strain, defined as e = {a^uh — ^sub)/(Q.sub + 
^sub), is negative at 100 K and increases upon heating. 



TABLE I. Atomic coordinates of a-(CuCl)LaNb207 at room 
temperature: a = 7.76290(8) A, 6 = 7.76197(7) A, c = 
11.73390(6) A, space group Pbam, and Ri/Rp = 0.019/0.063. 
Atomic displacement parameters Uiso for oxygen atoms were 
constrained. 
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FIG. 4. (Color online) Temperature dependence of 

the subcell lattice parameters and the subcell volume for 
(CuCl)LaNb207. Error bars do not exceed the size of the 
symbols. 



The strain amounts to zero at ~ 260 K and reaches 
its maximum value at 500 K. There is no indication 
of a structural transformation in the 100 — 500 K tem- 
perature range, thus the RT a-(CuCl)LaNb207 crystal 
structure is a reliable model for the investigation of the 
low-temperature properties. At 500 K, the orthorhom- 
bic strain starts to decrease, because of the a — >■ /3 
phase transition between the two orthorhombic poly- 
morphs. The strain vanishes at 640 K indicating a /3 — )■ 7 
orthorhombic-to-tetragonal phase transition. Although 
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FIG. 5. (Color online) Temperature evolution of the XRD 
patterns. The reflections are indexed in the 2aaub x 2asub x Csub 
supercell, only the superstructure reflections are labeled. The 
reflections with even h and odd k are attributed to the fci = 
[0, |,0] propagation vector, while the reflections with both h 
and k odd can be assigned to the fcg = [|,0,0] propagation 
vector. The patterns are offset for clarity. The maximum 
intensity for the main reflections is about 50000 counts. 



the two phase transitions are clearly visible in the SXRD 
data, we failed to observe them with differential scanning 
calorimetry (DSC), presumably, due to the low entropy 
change associated with these transitionsJ-^BJ The phase 
transitions are accompanied by a change in the intensities 
of the superstructure reflections (Fig. Isl). Above 640 K, 
the SXRD pattern of the 7-(CuCl)LaNb207 phase was 
indexed in a simple Ogub x cisub x Csub tetragonal unit 
cell. No reflection conditions were identified, suggesting 
the PA/mmm space group for the 7-phase. The SXRD 
pattern of the /3-(CuCl)LaNb207 phase shows additional 
reflections corresponding to the ki = [0, ij,0] propaga- 
tion vector. The resulting Osub x Sosub x Csub unit cell 
is orthorhombic, and the space symmetry is Pbmm {Okl, 
k — 2n reflection conditions). ESI The second set of su- 
perstructure reflections corresponds to the fc2 = [5,0,0] 
propagation vector^ and appears below 500 K at the 
/?—>■« phase transition. 

The continuous change in the unit cell volume (Fig. [4]) 
suggests that both phase transitions are of the second 
order. The PA/mmm — > Pbmm — > Pbam symmetry 
evolution is consistent with this assumption (see Sec. IV). 
The octahedral tilting distortion is completely suppressed 
in /3-(CuCl)LaNb207, and the superstructure is solely 
caused by the atomic displacements in the (CuCl) layer. 
We find ordered, " in-phase" shifts of the CI atoms along 
the a direction (see the middle panel of Fig. [6]) . The Cu 
atom can be placed into the special 2d position (^,0, ^), 
but the shift of Cu to a split, four- fold position reduced 
the refinement residuals and the atomic displacement pa- 
rameter of Cu. Therefore, the split position of Cu should 
be preferred. The refined atomic positions are listed 
in Table [nj The /3-(CuCl)LaNb207 structure essentially 



TABLE II. Atomic coordinates of /3-(CuCl)LaNb207 
at 500 K: a = 3.88666(3) A, 6 = 7.76396(5) A, 
c — 11.74441(7) A, space group Pbmm, and Ri/Rp = 
0.024/0.062. Atomic displacement parameters Uiso for oxy- 
gen atoms were constrained. 



Atom 


Position 


X 


y 


z 


a" 


4i 


0.5482(9) 





1 

2 


Cl 


2/ 


0.1145(9) 


3 

4 


1 

2 


La 


2e 


0.0001(3) 


3 

4 





Nb 


4/1 


1 
2 





0.80757(4) 


0(1) 


4ff 








0.8515(7) 


0(2) 


4fc 


0.508(2) 


1 
4 


0.8343(7) 


0(3) 


4/1 


1 
2 


1 

2 


0.6567(3) 


0(4) 


2c 


1 
2 









t/i. 



1.13(6) 
3.7(1) 
0.81(1) 
0.54(1) 
1.29(6) 
1.29(6) 
1.29(6) 
1.29(6) 



^ The f/iso values arc given in 10 ^ 
^ Occupancy factor 9=5 



TABLE III. Atomic coordinates of 7-(CuCl)LaNb207 at 
660 K: a = 3.88852(2) A, c = 11.73778(7) A, space group 
P4:/mm.m., and Ri/Rp = 0.027/0.066. Atomic displacement 
parameters Uiso for oxygen atoms were constrained. 



Atom 


Position 


X 


y 


z 


C/i.c" 


C,'^ 


4o 


0.5662(8) 


1 

2 


1 

2 


1.24(7) 


C^ 


4m 


0.096(2) 





1 
2 


4.7(2) 


La 


la 











1.08(2) 


Nb 


2h 


1 
2 


1 

2 


0.19305(4) 


0.69(1) 


0(1) 


Ai 





1 
2 


0.1582(2) 


1.88(6) 


0(2) 


2h 


1 
2 


1 

2 


0.3436(3) 


1.88(6) 


0(3) 


\c 


1 
2 


1 
2 





1.88(6) 



^ The Uiso values are given in 10^^ A^ 
^ Occupancy factor 9 = j 
'^ Occupancy factor 9 = j 



matches our DFT-based modeP^with respect to the unit 
cell and the crystal symmetry. The only difference is 
the partial disorder due to the split position of the Cu 
atom. In this position, each Cu atom has a short Cu~Cl 
bond of 2.34 A and a longer bond of 2.57 A. Similar to 
the Q;-(CuCl)LaNb207 structure, these bonds show the 
frans-configuration with a Cl-Cu-Cl angle close to 180° . 

In the tetragonal 7-(CuCl)LaNb207 structure, neither 
Cu nor Cl atoms lie on the four-fold axes. Our structure 
refinement (Table III and the right panel of Fig. [6| shows 
that both Cu and Cl partially occupy the split four- 
fold positions with a shortest Cu-Cl distance of 2.30 A. 
This implies that copper prefers to keep short bonds to 
the chlorine atoms, but the spatial arrangement of these 
short bonds is now fully disordered. 
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FIG. 6. (Color online) The structure of (CuCl)LaNb207 at different temperatures. Left panel: an overall view with the tilted 
NbOe octahedra and the fully ordered [CuCl] layer below 500 K (a-modification). Middle panel: the [CuCl] layer with the 
disordered Cu atoms (between 500 K and 640 K, /3-modification) . Right panel: the disordered structure of the [CuCl] layer 
(above 640 K, 7-modification). The arrows show the primary atomic displacements. Note the interconnected displacements of 
Cu and O in the left panel. 



C. Electron diffraction and transmission electron 
microscopy 

Electron diffraction patterns of a-(CuCl)LaNb207 
(Fig. [t]) were readily indexed in the pseudo-tetragonal 
unit cell with a ~ 6 ~ 7.75 A and c ~ 11.8 A, as proposed 
by the SXRD data. These patterns clearly evidence the 
superstructure in the ab plane. The [001] ED pattern in 
Fig. [Tk has apparent tetragonal symmetry. However, a 
closer inspection of different regions of the crystals with 
a smaller selected area aperture revealed the areas where 
the intensity distribution contradicts the four-fold sym- 
metry (Fig. [7b). The pseudotetragonal [001] ED pattern 
(Fig. I7k) is then a sum of the ED patterns from 90° rota- 
tional twin domains with orthorhombic symmetry. Due 
to the small difference between the a and b lattice pa- 
rameters and due to the twinning, the [010] and [100] 
ED patterns could not be distinguished (Fig. ITb). The 
reflection conditions h = 2n for hOl and k ^ 2n for Okl, 
derived from the ED patterns, are in agreement with the 
Pbam space group. In the [001] ED pattern, the forbid- 
den OfcO, k ^ 2n reflections appear due to the multiple 
diffraction. 

We found that the ordering in the (CuCl)LaNb2 07 
structure is very sensitive to electron-beam irradiation. 
The superstructure reflections with h + k ^ 2n disappear 
rapidly and irreversibly under an intense electron beam, 
while the reflections with h + k — 2n remain. HRTEM 
images demonstrate that these changes are related to the 
amorphization of the (CuCl) layers under the electron 
beam, whereas the LaNb2 07 blocks are less affected. The 
evolution of the structure under the electron-beam irra- 
diation is surprisingly different from the effect of heat- 
ing (Sec. IIIB). The heating first removes the tilts of 
the NbOg octahedra, while the (partial) ordering in the 
(CuCl) layers survives up to 640 K. In contrast, the elec- 
tron beam "melts" the (CuCl) layers and weakly changes 
the perovskite blocks. Weak diffuse intensity lines along 
c* are visible on the [010]/[100] and [110] ED patterns 
(Fig. [7b and d) and can be attributed to a local irradia- 
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FIG. 7. Electron diffraction patterns of a-(CuCl)LaNb207. 
The superstructure reflections are present in the patterns a, b, 
and d. The pattern a shows the pseudotetragonal symmetry 
due to the 90° rotational twin domains, while the pattern b is 
taken with a smaller aperture and points to the orthorhombic 
symmetry. White boxes denote several positions of the miss- 
ing reflections. The patterns c and d mainly show the subcell 
reflections. 



tion damage of the (CuCl) layers. 

Remarkably, the HAADF STEM mode does not intro- 
duce such a severe beam damage as the HRTEM imaging. 
On the HAADF STEM image, the intensity is related 
to the average atomic number of an atomic column as 
/ « Z" (n = 1 - 2). The [100] HAADF STEM image of 
(CuCl)LaNb207 (Fig. [8| reveals a defectless sequence of 
the cationic layers along the c axis. However, the oxygen 
and chlorine atoms have too low atomic numbers to be 
observed in this imaging mode. The oxygen and chlorine 




FIG. 8. [100] HAADF STEM image of a-(CuCl)LaNb207. 
The cation positions are overlaid. 






tt»»t«w M wH»*> 



r.i.v.i:i:i:i"'.r.:rr.ti:Kiiir.v.iunana:Tx:r 




FIG. 10. [001] HRTEM image of a-(CuCl)LaNb207. The 
Fourier transformed image is shown at the top left corner. 
The simulated image (/ — —38 nm and t — 3.5 nm) is shown 
as an inset and outlined with a white rectangle. Arrowheads 
mark the rows of brighter dots. 
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FIG. 9. [110] HRTEM image of a-(CuCl)LaNb2 07. The 
simulated image (/ = 31 nm and t = 3.3 nm) is shown as an 
inset and outlined with a white rectangle. 



Fourier transformed image, which shows the spots corre- 
sponding to the doubling of the lattice parameters in the 
ab plane. The simulated HRTEM image (/ = —38 nm, 
t = 3.5 nm) calculated with the Pbam model lib is in 
agreement with the experimental data. 



D. Structure relaxation and electric field gradients 



displacements can be revealed by the HRTEM technique. 
In order to reduce the beam damage, the HRTEM images 
were collected under the minimum possible beam inten- 
sity and keeping the exposure time as short as possible. 
On the [110] HRTEM image (Fig. [9]), the rows of brighter 
dots correspond to the empty spacings in the layers of 
apical oxygen atoms of the NbOe octahedra. The dots 
within the brighter rows are arranged into alternating 
pairs with different brightness due to the displacements 
of the apical oxygen atoms upon the octahedral tilt. The 
theoretical HRTEM image (defocus / = 31 nm, sample 
thickness t = 3.3 nm) calculated with the Pbam model 
lib is in agreement with the contrast of the experimen- 
tal image. Note, however, that the displacements in the 
(CuCl) layers can not be observed in this structure pro- 
jection. 

The atomic displacements in the (CuCl) layer are vis- 
ible when the structure is viewed along the [001] axis 
(Fig. 10 1. Under these imaging conditions, the brighter 



dots correspond to the projections of the LaCl columns. 
The superstructure is clearly visible on this image as an 
alternation of brighter and less bright rows (some of the 
brighter rows are marked with arrows in Fig. 10 1. The 



presence of this superstructure is also confirmed by the 



The experimental information on the unit cell and 
the crystal symmetry imposes helpful constraints for the 
computational analysis. The calculations, on the other 
hand, can evaluate the relative stability of different struc- 
tural models. If all the atomic coordinates are fixed to 
the experimental positions, model lib shows the low- 
est energy. Models Ha and HI lie higher in energy by 
0.53 eV/f.u. and 0.81 eV/f.u., respectively. As the 
atomic positions are relaxed, models Ha and lib con- 
verge to the same solution that closely resembles model 
lib with short Cu-Cl bonds in the trans-configuration. 
The relaxation of model HI yields a similar structure 
with slightly higher energy, because the Cu atoms are 
constrained to the two-fold positions. Overall, the com- 
parison of total energies and the ab initio structure relax- 
ation strongly support model lib and discard the other 
two models of the a-(CuCl)LaNb207 structure. 

It is also instructive to compare the details of the 
experimental and relaxed geometries. Table |IV| shows 
the interatomic distances and the Cu-Cl-Cu angle, 
as found from the structure refinement and from the 
LSDA-|-[//GGA-t-l7 relaxations. The LSDA-^C/ re- 
laxation underestimates the Cu-Cl distances, while 
GGA+U shows better agreement with the experiment. 



TABLE IV. Interatomic distances (in A) and angles (in deg.) 
in Q:-(CuCl)LaNb207, as found from the structure refinement 
and from the LSDA+L*' or GGA+U ab initio relaxations. 



Refinement 



LSDA+C/ 



GGA+U 



Cu-Cl 

Cu-0(3) 

Cu-Cl-Cu 

Nb-O(l) 

Nb-0(2) 

Nb-0(3) 
Nb-0(4) 
Nb-0(5) 



2.41(1) 
2.43(1) 
2 X 1.85(1) 
107.1(2) 
1.96(1) 
1.97(2) 
2.02(2) 
1.78(1) 
2.26(1) 
2.02(1) 



2.32 
2.32 
2 X 1.85 
113.7 
1.97 
1.99 
2.00 
1.79 
2.26 
2.03 



2.38 
2.38 
2 X 1.85 
109.7 
1.93 
2.00 
2.00 
1.79 
2.29 
2.04 



Thus, GGA+U is the best available computational ap- 
proach to model the crystal structure of (CuCl)LaNb207. 
Our refined interatomic distances are in remarkable 
agreement with the previous structure refinement based 
on the neutron diffraction data.'i^ The neutron experi- 
ment did not show any superstructure reflections, and 
the ordered arrangement of the CI atoms could not be 
found. However, the CI atom was moved to the four- 
fold position with an occupancy of j, yielding a Cu-Cl 
distance of 2.40 A (2.41 - 2.43 A in our refinement). Ac- 
cording to Ref. lHj the Cu-0 distance is 1.84 A; this 
also matches our findings (1.85 A). The reference to the 
neutron data ensures the remarkable accuracy of the ob- 
tained structural information, despite the huge difference 
in the scattering power of oxygen and the heavy (La, Nb) 
atoms. 

Structural distortions can be traced by EFGs that have 
been determined in a NMR experiment .1^^ Table \V\ lists 
EFGs calculated for three different models of the RT 
structure: Vzz is the leading principal component of the 
tensor, while r] = (Vyy — Vxx)/Vzz is the asymmetry. The 
results for the Cu and CI sites are similar to our previous 
report: the z axis of the EFG tensor lies in the ab plane 
for Cu and aligns along c for Cl.^^ Note however that 
the model Ila leads to a pronounced asymmetry of the 
Cu EFG (ry = 0.43) due to the very irregular local envi- 
ronment of the Cu atom. The model lib shows a lower 
asymmetry (77 = 0.04) in agreement with the experimen- 
tal result {rj — 0.10). Finally, the model III yields fairly 
different EFGs (well above the experimental resolution) 
for the two inequivalent Cu sites and is readily discarded. 

The tilts of the NbOg octahedra allow to reproduce 
the sizable asymmetry of the EFG at the La site. The 
calculated asymmetry for the model lib {rj — 0.57) is in 
agreement with the experimental value (77 — 0.70). In 
contrast, the model Ila yields a very low asymmetry and 
does not fit to the experimental findings. Compared to 
the experimental data, the V^z value remains rather high 
in all the three models. This discrepancy can be caused 



TABLE V. Calculated EEC's for three models of the a- 
(CuCl)LaNb207 structure. Vzz is the principal component of 
the tensor (in 10^^ V/m^), 77 is the tensor asymmetry. Model 
III has two inequivalent Cu positions. The experimental val- 
ues are taken from Ref. T^ and do not allow to determine the 
sign of Vzz- The structural models are shown in Fig. [2] 



Cu 



CI 



La 



Vz. 



Vz, 



V 



Vz. 



V 



Model Ila 

Model lib 

Model III 

Experiment 



-12.9 0.43 

-11.4 0.04 

'.2/-12.2 0.16/0.22 

11.6 0.10 



-16.6 0.73 -17.3 0.08 

-18.3 0.50 -14.2 0.57 

-18.0 0.44 -21.2 0.26 

-14.2 0.56 6.5 0.70 



by 4/ orbitals of La or Ad orbitals of Nb that lie close 
to the Fermi level and are subject to correlation effects. 
Our present approach accounts for the correlations in the 
Cu 3d shell only and treats all the other orbitals on the 
LDA level. To achieve better agreement with the exper- 
iment, more sophisticated computational techniques are 
necessary. Apart from the deviation for Vzz, the calcu- 
lated EEC's strongly support the model lib and clearly 
disfavor the two other models. 

The ab initio structure relaxation also provides insight 
into the interplay of the fci and fc2 superstructures in a- 
(GuCl)LaNb207. To reveal the influence of the tilting 
distortion on the structure of the (CuCl) layers, we per- 
formed a relaxation with fixed positions of La, Nb, and O. 
We used the regular structure of the (LaNb207) blocks, 
as found from the refinement at 660 K. Taking the Cu 
position as {j+x,^+y, ^), we find a negligible displace- 
ment of the Cu atoms {x = -0.0001 and y = 0.0001) for 
the relaxed structure with the fixed (LaNb207) block, 
compared to a; = —0.0157 and y — 0.0058 for the ex- 
perimental crystal structure. Thus, the tilting distortion 
tends to shift the Cu atoms from their "parent" (|, ^, ^) 
position. In the Q;-(CuCl)LaNb207 structure, the octahe- 
dra tilt around the b axis, hence the axial 0(3) atoms (see 
Fig. fTl) displace along the a axis. The Cu atoms follow 
these displacements and also shift along a (left panel of 
Fig. [6]) . The interconnected displacements of the Cu and 
0(3) atoms underlie the coupling between the fci and ^2 
superstructures in (CuCl)LaNb207. This coupling can 
be understood as the tendency to keep the regular four- 
fold coordination of Cu with Cu-0 bonds perpendicular 
to the Cu-Cl bonds. The preference to the regular co- 
ordination is also evidenced by the trans-arrangement 
of the Cu-Cl bonds and rules out alternative structural 
models that would suggest a highly irregular local envi- 
ronment of the Cu atom. 



IV. DISCUSSION 

In the preceding section, we have derived the a- 
(CuCl)LaNb207 structure by considering different order- 
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ing patterns in the (CuCl) layer (Fig. [2]) and selected 
the model lib, based on the comparison of the refine- 
ment residuals, the total energies, and the EFGs. Alter- 
natively, one can consider the a-(CuCl)LaNb207 struc- 
ture as the derivative of the high-temperature tetrag- 
onal 7-(CuCl)LaNb207 phase. Assuming second-order 
phase transitions, the structural distortion in the low- 
temperature polymorph should follow one of the irre- 
ducible representations for the space group of the high- 
temperature polymorph. 

The Pbmm /3-structure can be derived from the high- 
temperature PA/mmm 7-structure by the X^ irre- 
ducible representation [here X denotes the (0, |,0) point 
in the reciprocal space and implies a doubling of the b pa- 
rameter]. This irreducible representation decouples the 
four CI atoms within one position and splits the four 
equivalent Cu atoms into two pairs. The leading dis- 
tortion is the "in-phase" shift of the CI atoms. The 
Cu atoms are disordered over two equivalent positions 
around (5,0, ^). 

The /3 — >■ a transition follows the Z^ irreducible repre- 
sentation of the Pbmm space group, where Z is (i, 0, 0) 
and corresponds to a doubling of the a parameter [Z is 
(0,0, |) in the standard Pmma setting]. The Cu atoms 
shift along a in the " out-of-phase" manner (the middle 
panel of Fig. l6|. The 0(3) atoms also shift along a, 
leading to the a^b~c'^ octahedral tilting pattern. The 
structural changes at 500 K and 640 K can be identi- 
fied as atomic displacements order-disorder phase transi- 
tions. The evolution of the crystal symmetry is in agree- 
ment with the second-order nature of these transitions. 
Note that the intermediate Pbmm structure allows for 
the "in-phase" shifts of the CI atoms only, thus impos- 
ing the similar arrangement of the CI atoms in the a- 
structure. The alternative structural model by Yoshida 
et ali^ (model Ila) would require a dramatic rearrange- 
ment of the (CuCl) layers, which is unlikely. 

Since the experimentally determined structure of the 
(CuCl) layer in the model lib is rather similar to 
our previous DFT-based model,'^^ certain conclusions 
on the electronic structure and magnetic properties 
can be derived. The regular (tetragonal) structure of 
(CuCl)LaNb2 07 leads to the orbital degeneracy of Cu. 
Atomic displacements in the (CuCl) layers are primarily 
caused by the tendency to lift the orbital degeneracy. The 
resulting coordination of the Cu atoms is the CUO2CI2 
plaquette with the magnetic (half-filled) dx2-y2 orbital 
lying in the plane of the plaquette. The respective bands 
are well separated from other Cu 3d states and ensure 
the one-orbital scenario, typical for cuprates. 

The lobes of the magnetic dx^-y'^ orbital point toward 
CI and O atoms. The strong cr-overlap between Cu 3d 
and CI 3p orbitals leads to a rather unusual superex- 
change scenario with strong coupling between fourth 
neighbors (Cu-Cu distance of about 8.7 A) in agreement 
with the inelastic neutron scattering results.™ A quan- 
titative evaluation of this unusual long-range coupling 
depending on the CI position has been given in Ref. Il7l 



The strong hybridization of the magnetic Cu 'id^^-y^ 
orbital with the oxygen orbitals leads to a somewhat 
three-dimensional nature of the system. In Ref. [T71 we 
have shown that this hybridization along with the low- 
lying Ad states of Nb results in a sizable interlayer cou- 
pling of about 15 K, compared to 50 — 60 K for the cou- 
plings in the ab plane. The strong hyperfine coupling at 
the Nb site^^ can be taken as an experimental signature 
of the strong Cu-O-Nb hybridization. 

Although the reference to the DFT-based model clari- 
fies the electronic structure of the compound, we should 
preclude the reader from the straight-forward transfer 
of the microscopic magnetic model. Superexchange cou- 
plings are highly sensitive to fine details of the crystal 
structure. Therefore, sizable changes in the exchange 
couplings are expected. According to Table llVl the 



LSDA-|-[/ structure relaxation leads to a Cu-Cl-Cu angle 
of 113.8°, while the experimental value is about 107.1°. 
The angles at the ligand atoms are known to have a dra- 
matic influence on the magnitude and even on the sign 
of the superexchange couplings. Following simple micro- 
scopic arguments in the spirit of Goodenough-Kanamori 
rules,'^we suggest that the AFM nearest-neighbor inter- 
action along the b direction should be reduced compared 
to our previous calculationgi^ or even changed towards a 
FM interaction. Since this interaction is the leading cou- 
pling in the original spin model, substantial changes in 
the resulting magnetic ground state should be expected. 
The shifts of the Cu atoms result in a large number of in- 
equivalent Cu-Cu distances and, consequently, in a large 
number of different exchange interactions. The detailed 
analysis of the spin model is a challenging task involving 
extensive band structure calculations, model simulations, 
and quantitative comparison to the experimental data. 
Such a study is presently underway.'221 

Using the combination of different techniques, the 
present study gives a detailed insight into the chemi- 
cal modifications of (CuCl)LaNb207. We show that the 
atomic displacements in this crystal structure are inter- 
connected, and a tilting distortion of the perovskite block 
induces atomic displacements in the (CuCl) layers (see 
the left panel of Fig. [6]). This explains why the substi- 
tutions in the perovskite block (e.g., Ta for Nb) have 
pronounced effects on the magnetism,E2l despite the lead- 
ing interactions are expected within the (CuCl) layers for 
the whole compound family. The interplay of the atomic 
displacements also provides a key to the puzzling behav- 
ior of (CuBr)Sr2Nb30io and related compounds where 
tiny changes in the pero vskite block destroy the peculiar 
I -magnetization plateau.ESEll 

Based on the temperature evolution of the 
(CuCl)LaNb207 structure, we are able to propose 
a general mechanism for structural distortions in the 
family of Cu-based quantum magnets as well as in 
isostructural (MCl)LaNb207 compounds (M = Cr, Mn, 
Fe, Co).'2SH2Zl The regular (tetragonal) structure of the 
[CuX] (X = CI or Br) layer is highly unfavorable due to 
the orbital degeneracy of Cu (Ref. IT7|) and unsuitable 
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Cu-X distances. The strain is released by the displace- 
ments of the Cu and X atoms leading to shorter and 
longer Cu-X bonds. Below certain temperature [640 K 
in (CuCl)LaNb207], such shifts become cooperative, but 
the Cu atoms remain disordered due to their bonding to 
the oxygen atoms (middle panel of Fig. [6]) . The tilting 
distortion allows to shift the Cu atoms and results in 
their ordered arrangement (left panel of Fig. ^. Thus, 
the low-temperature structure is fully ordered and can 
be reliably determined using conventional diffraction 
techniques. The interplay of the atomic displacements 
plays a crucial role in this compound family. The 
specific tilting distortion sets an ordered configuration 
of the Cu atoms and has a strong effect on the magnetic 
properties. Further studies of related materials should 
challenge this mechanism and disclose the influence 
of different atoms (CI vs. Br and Nb vs. Ta) on the 
transition temperatures and distortion patterns. 

In summary, we have determined the room- 
temperature crystal structure of (CuCl)LaNb2 07 and ob- 
served two structural order-disorder phase transitions in 
this compound. The room-temperature structure com- 
bines the tilting distortion of the (LaNb207) perovskite 
blocks with ordered displacements of the Cu and CI 
atoms within the (CuCl) layers. The tilting distor- 
tion disappears above 500 K and leads to the partially 
disordered arrangement of the Cu atoms while a fur- 
ther transition at 640 K destroys the ordering in the 
(CuCl) layers completely. The room-temperature struc- 
ture shows a clear preference of Cu to the four-fold 
square-like coordination with two O and two CI neigh- 
bors in a CUO2CI2 plaquette. The plaquettes share cor- 
ners and form chains running along the b direction. This 
type of the structure organization is basically consistent 



with our previous computation-based structural model 
and resolves the controversy between several compet- 
ing proposals. The obtained structural information is 
the necessary and long-sought input for understanding 
(CuCl)LaNb207 and related compounds with unusual 
magnetic properties. 

Note added: after the submission of the original 
manuscript, the results of Ref.lT51were updated.l^The re- 
ported space group and the DFT-based structural model 
are in agreement with our experimental structure solu- 
tion, although the interatomic distances and angle s are 
slightly different, similar to Table IV. Ren and Chen^iS138| 
accidentally refer their model to the earlier proposal by 
Yoshida et a/., (Ref. 1I61) despite the corner-sharing ar- 
rangement of the CUO2CI2 plaquettes essentially matches 
our original model ('Ref. [T7|) t hat was also based on DFT. 
Additionally, Tassel et al^ refined the structure of a- 
(CuCl)LaNb207 using neutron diffraction data and also 
arrived to the orthorhombic Pbam symmetry. The re- 
sulting atomic positions are in close agreement with our 
results. 
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